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Impurity Site NMR Relaxation in Unconventional Superconductors 
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Impurity nuclear spin relaxation is studied theoretically. A single impurity generates a bound 
state localized around the impurity atom in unconventional superconductors. With increasing 
impurity potential, the relaxation rate T^^ is reduced by the impurity potential. However, it 
has a peak at low temperatures due to the impurity bound state. The peak disappears at non- 
impurity sites. The impurity site NMR measurement detecting a local electronic structure just 
on the impurity atom is very useful for identifying the unconventional pairing states. 
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The study of unconventional superconductors has be- 
come one of the most attractive issues in recent con- 
densed matter physics. These materials include heavy 
fermion compounds, high- Tc cuprates, organic supercon- 
ductors and the recently discovered Sr2Ru04. A great 
effort has been made to detect the unconventional su- 
perconductivity by several probes: specific heat mea- 
surement, tunneling spectroscopy, nuclear magnetic res- 
onance (NMR), TT-junction and so on. NMR experiments 
were important in the early stage of studying the uncon- 
ventional superconductivity. It is well known that T^^ 
in s-wave superconductors has a Hebel-Slichter peak just 
below Tc, while it does not exhibit this feature in uncon- 
ventional superconductors. On the other hand, it was 
also reported that the peak can be suppressed by the 
strong damping effect even for a s-wave state. tP 

To clarify this point, Ishida et al. performed a 
Cu site (bulk site) NMR experiment with Zn-dopcd 
YBaaCusOy-s (YBC0)J3) They found that Tf ^ is pro- 
portional to T at low temperatures. Non-magnetic impu- 
rities do not break the s-wave superconductivity.EP How- 
ever, they break unconventional ^.upierconductivity and 
yield low-energy impurity statesD'aEP The observed T 
linear dependence of T^^ can therefore be understood 
as evidence of extended low-energy impurity states and 
supports an unconventional pairing state for YBCO. In 
this case, the effect of many impurities was important. 

A single impurity can also break the unconventional 
superconduatwity locally and generate a bound state 
around it.Q^Q'LP Recently, Pan et al. observed such a 
low-energy state around a Zn atom in Bi2Sr2CaCu208+5 
using a low-temperature scanning tunneling microscope 
(STM), showing a fourfold symmetrical structure whi 



is consistent with a dx2-y2-wave pairing state£3'La 
Nishida et al. also observed the fourfold structure around 
a columnar defect by a STmIb^ At present, STM is 
the only instrument which can reveal the local electronic 
structure around a single impurity atom. 

Usually, NMR is used to probe bulk electronic struc- 
ture. For example, Cu site NMR measurements for high- 
Tc cuprates detect local electronic structure at the Cu 



atom. Recently, impurity sitaJ^MR experiments were 
performed by several groups.llalliP However, there is, so 
far, no theory for the impurity site NMR. The purposes 
of this paper are (1) to address the theory for the im- 
purity site NMR and (2) to propose a new NMR experi- 
ment which probes the local quasiparticle states around 
a single impurity atom. In this study we concentrate our 
attentions on a single impurity problem, since we can 
solve it exactly. For unconventional superconductivity, 
we show that the impurity site NMR T^^ has a peak, 
while the peak disappears at the bulk site (non-impurity 
site). This peak at the impurity site is a common feature 
of all types of unconventional superconductors. The im- 
purity site NMR can reveal the local quasiparticle struc- 
ture as well, and provides us with a new experimental 
method for identifying unconventional superconductiv- 
ity. Recentlv_MMR T^^ in a vortex state was studied 
theoreticallyJlaM They proposed that the NMR can be 
used as a site-selective probe by controlling the resonance 
frequency. Thus, NMR is gradually attracting attention 
as a local probe. 

In the unconventional pairing case, the superconduct- 
ing order parameter is suppressed around the impurity. 
However, we present our theory assuming a uniform or- 
der parameter, since we can capture the essential physics 
clearly without any complex analysis. Actually, we have 
performed numerical self-consistent calculations to in- 
clude the spatial dependence of the order parameter. We 
have confirmed that the NMR T['^ has no qualitative dif- 
ference between in the uniform and non-uniform order 
parameter cases. The most important point in the im- 
purity site NMR is the existence of the impurity bound 
state, which does not depend on details of unconven- 
tional pairing states. 

For the unconventional pairings, we focus on px ± ipy- 
wave and dx2_y2-wave pairing states. The p^; ± ipy-wave 
is the simplest and most essential, sjojoraetry suggested 
for the Sr2Ru04 superconductor .Il3ll§ll3) This state is a 
spin—tpplet state and breaks the time-reversal symme- 
try.E3'ciP For simplicity, we assume that the supercon- 
ductor is basically two-dimensional and has a cylindrical 
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Fermi surface. The Mataubara Green function in a 2 x 2 
matrix form is given byEP 

G{iujm,r,r') = Go(it^m,r,r') + Go(icj„i, r, 0)C/of3 

1 - Go(lWm,0,0)f/oT3 

where a single impurity is located at the origin of the co- 
ordinate, [/o represents the strength of the short-range 
impurity potential. u)m is the fermion Matsubara fre- 
quency and Ti {i — 1, 2, 3) represents the Pauli matrix in 
a charge space. For simplicity we use h = 1 and ks — 1 
units throughout this letter. Go is the non-perturbed 
{Uo = 0) Green function given by 



Go{iuj„i,r,r') = 



k 



ik-{r-r')] 



V3 



(2) 

Here, is equal to A^fi and expresses the momentum- 
dependent order parameter. For the px ± ipj,-wave 
(d,2_y2-wave), Afe = Ap(r)e"''' [A^ = Ad(T) cos(2^fc)]. 
Here, 6k is the angle of the Fermi wave vector measured 
from the fc^ axis. We solve the gap equation in the bulk 
region and use the temperature-dependent order param- 
eter Ap(T) (Ad(T)) which is real and positive for the 
Px ± ipy-wave ((i^2_j,2-wave). is the kinetic energy 
and Vl represents the volume of the system. NMR relax- 
ation occurs via the hyperfine interaction between the 
nuclear spin and the conduction electrons. I vlMR Tf ^ at 
the impurity site is given in an explicit forrr 



12) 



Tfi =2^( — )2(7,7„)2h/, 
,aii{E)a22{-E) 



W ^ AE 



ai2{E)a2i{-E) 



1-f cosh(£;/T) ' 

where W is proportional to the nuclear spin flip tran- 
sition probability, je and 7„ are the gyromagnetic ra- 
tios for the electron and nucleon, respectively. (E) = 
— Im[Gij(iWm ^ E + i6, 0, 0)]/7r, where the subscript of 
ttij represents the matrix element of the Green function. 
Here, we have omitted the uq dependence in for sim- 
plicity, where uq = ttNqUo (Nq is the density of states 
per volume at the Fermi energy). (5 is a positive, small 
number which expresses a finite level broadening. 

In the s-wave case, NMR Tf ^ has the Hebel-Slichter 
peak just below Tc at the non-impurity site [uq — 0). 
At the impurity site [uq ^ 0), aij scales as aij{E,uo) — 
aij{E, 0)/(l + Uq). The local density of states at the im- 
purity is then simply reduced, keeping the same energy 
dependence. This is consistent with Anderson's theo- 
rem.cP In the same manner, W scales as 



W{uo,T)^W{Q,T)/{l + uiy^ 



(4) 



Therefore, 



is simply reduced at the impurity site, 



while the temperature dependence does not change be- 
tween the impurity site and the bulk site for the s-wave, 
as shown in Fig. 1(a). 

For the p^, ±ipy-wave, we show the results in Fig. 1(b). 
In this case, the temperature dependence at the impu- 
rity site is quite different from that of the bulk site. In 
the unconventional pairing case, both ai2 and a2i are 
zero, since pair electrons cannot possess the same posi- 
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Fig. 1. Temperature dependence of W with various uq with a 
fixed damping rate S/Tc = 0.1. The number represents a value of 
Uq. (a) and (b) are for the s-wave and pxiipy-wave, respectively. 
The peak of Tj~^ in (b) becomes higher with the decrease of 5/Tc 
[see eq. (|^)]. The value of 5 /Tc docs not change the temperature 
dependence of ^ . 



tion. Therefore the coherence factor vanishes for uncon- 
ventional pairings, resulting in the absence of the Hebel 
Slichter peak. W is then expressed by the following sim- 
ple form: 



W = N, 



dE- 



where Ni„,p{E) 



l + cosh{E/T)' 

an{E)/No = a22i-E)/No is the di- 
mensionless local density of states at the impurity atom. 
For the Px ± ipj,-wave, a bound state appears around the 
single impurity. The energy position of the bound state 
is given by Eb = -sgn{uo)Ap{T) / ^/TTv^-^ In fact, 
an impucity effect is observed in Sr2Ru04 as a reduction 
of its TcE5> iVinip in eq. (||) for the px ± ip^-wave has the 
following analytic form: 

^'_M^S(E-E,). (6) 

(1 + M^)2 



where Np{E) = \E\/Je^ - A2(r) is the dimensionless 
Px ± ipj,-wave bulk density of states. 6 and S in eq. 
(6) are the Heaviside and 5-functions, respectively. The 
first and second terms in eq. (6) are the continuum and 
bound states, respectively. Introducing uq, we can see in 
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Fig. 2. (a) Dimensionless local density of states for the ■px itipy- 
wave state. TVp is multiplied by (1 + = N{^p{oo) for 

convenience. Parameters are chosen as ttg = 1 and 5/Ap = 0.01. 
(b) Temperature dependence of Wb- Parameters are chosen as 
uo = l and 5/Tc = 0.1. 



Fig. 2(a) that the density of states from the quasipar- 
ticle continuum is decreased, while that of the impurity 
bound state increases. As in the s-wave case, W for the 
Px ± ipy-wave state is reduced with the increase of uq. 
However, as shown in Fig. 1(b), it drasticaUy changes at 
low temperatures due to the impurity bound state. The 
bound state contribution to W behaves as 
Ap{T)uQi^ 



(1 



cosh{EB/T)]- 



(7) 



where 6{E - Eb) = {S/tt)[{E - EbY + 5^]"^ has been 
used. Figure 2(b) clearly shows that Wb has a peak be- 
low Tc- Notice that the peak in W in Fig. 1(b) is closely 
related to the local impurity bound state. It is different 
from the Hebel-Slichter peak and can be distinguished, 
since T^^ does not exhibit such a peak at the bulk site 
for unconventional superconductivity. The peak appears 
only at the impurity site. At low temperatures, impurity 
nuclear spin relaxation occurs via the impurity bound 
state at the impurity site. As uq is increased, the peak 
position of W shifts to a lower temperature region [see 
Fig. 1(b)], which reflects the energy shift of the bound 
state. For a large wq, Wb can survive if 5 is sufficiently 
small. 

Next we study the da.2_y2-wave case, which is the most 
favorable pairing symmetry for high-T,, cuprates. In this 
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Fig. 3. Temperature dependence of W with various uq for the 
(^2^2 _y2-wave with a fixed damping rate <5/Tc = 0.1. The number 
represents a value of uq. 



case, iVimp in eq. (||) takes the following form: 

1 Nd{E) 



ul [A{E) 



where A{E) and Nd{E) i 
nary parts of ((i? + i^)/-\/( 



\lu^\^-rNl{Ey 

e the real and 
-ii? + <5)2 + |AfcP) 



(8) 



imagi- 
re- 



fcl /FS' 

spectively. Here, (• ■ •)fs represents an average over the 
Fermi surface. The energy position of the bound state is 
determined by the zero of A{E) + I/wq and its effective 
broadening is Nd{E) which is the dimensionless d^'^-y'^- 
wave bulk density of states. Similar to the Px ± ipi/- 
wave case, the impurity site (uq 7^ 0) NMR Tj^^ for the 
d^2_j,2-wave has a peak due to the impurity bound state, 
as shown in Fig. 3. 

The important point in this letter is that a single im- 
purity generates a bound state in unconventional super- 
conductors, and that the NMR experiment can detect 
the local impurity bound state by the impurity site T^^ 
measurement. We emphasize that the peak in T^^ at 
low temperatures in the impurity site NMR is a com- 
mon feature exhibited by all types of unconventional su- 
perconductors. It does not depend on the details of un- 
conventional pairing states. We have demonstrated this 
by examining the Px ± ipj,-wave and dj.2_j^2-wave pair- 
ing states. The former is a fully gapped unconventional 
state, while the latter is a gapless one. 

Recently Ishida et al. carried out NMR and NQR mea- 
surements on CeRu2 with 1% Al (impurity) substitution 
for Ru.y They reported that the Tf ^ at the Al site 
was reduced to 10% of that of the Ru (bulk) site. The 
strength of the impurity scattering can be estimated as 
I Mo I = 1.5 by using eq. (Q). If the pairing is unconven- 
tional, we can expect the appearance of a peak below 
Tc for |mo| = l-5j as shown in Figs. 1(b) and 3. How- 
ever, the Al site Tj~^ does not exhibit such a peak. It is 
simply reduced, keeping the same exponential tempera- 
ture dependence as that of the Ru site. This feature is 
consistent with our theory of impurity site NMR for the 
s-wave. 

Very recently, Bobroff et al. succeeded in Li (impurity) 
site NMR measurement in YBCO with Li substituted for 
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Cu.M The Li concentrations of the samples were 0.85% 
and 1.86% per Cu02 layer. For YBCO, Li induces a mag- 
netic moment around it. Unfortunately, it is difficult to 
extract the contribution from the conduction electrons in 
Ti, since the induced moment also generates the relax- 
ation. However, this Li site experiment shows that the 
recent NMR experiment has enough sensitivity to detect 
Ti at the impurity site in unconventional superconduc- 
tors. The impurity concentration is about 1% for both 
CeRu2 and YBCO. These experiments indicate that the 
impurity site NMR is possible at least for such impurity 
doping. 

For the bulk site NMR Ti measurement, it is im- 
portant to observe whether the Hebel-Slichter peak is 
present or not. Examining the temperature dependence 
(exponential or power) at low temperatures is also im- 
portant. (I) In the absence of the Hebcl-Slichtcr peak, 
it is difficult to distinguish unconventional pairings from 
the s-wave for the bulk site measurement, since strong 
damping can suppress the coherence peak for the s-wave. 
(II) The temperature dependence of the bulk site T^^ is 
not useful for distinguishing fully gapped unconventional 
states from the s-wave. We cannot distinguish gapless 
unconventional states from the anisotropic s-wave, ei- 
ther. In contrast to the bulk site measurement, the im- 
purity site measurement can distinguish for both (I) and 
(II) cases. The feature of unconventional pairings ap- 
pears as a peak in T^^ only at the impurity site, while 
the peak disappears at the bulk site. On the other hand, 
for the s-wave, the temperature dependence of T^^ does 
not change for either impurity or bulk sites [see eq. (^) 
and Fig. 1(a)]. If we combine the two (bulk and impurity 
sites) NMR measurements, we can identify the uncon- 
ventional pairing without ambiguity. Thus, the impurity 
site NMR is a powerful experimental method for probing 
unconventional superconductivity. 

Our theory for a single impurity is applicable to a case 
with low impurity concentration. With the increase of 
the impurity concentration, extended impurity states are 
formed by the local impurity bound states, resulting in 
additional continuum energy levels below the energy gap. 
In this case, the impurity effect appears as a constant 
TiT at low temperatures for both the impurity and bulk 
sites. However, the impurity site measurement is ex- 
pected to exhibit the impurity effect more obviously than 
the bulk site, since the electronic states at the impurity 
site are strongly affected by the impurity atom. We note 
that the peak in T^^ of our theory for the single impu- 
rity (low impurity concentration) is much more drastic 
than the constant TiT, both of which are evidence of 
unconventional superconductivity. 

The local electronic structure around the impurity 
atom reflects the superconducting pairing symmetry. So 
far, STM is the only probe of the local impurity state. 
However, STM experiments are strongly affected by the 
surface roughness. We note that the impurity site NMR 
can be used as a local probe regardless of annoying sur- 
face conditions. 

In conclusion, we have addressed the theory for the 
impurity site NMR T^^ in unconventional superconduc- 
tors. We have concentrated our attention on a single 



impurity effect. There are two characteristic points of 
the impurity site NMR. The first one is that the relax- 
ation rate Tj~^ is reduced by the impurity potential. The 
second is the peak in T^^. It appears as direct result 
of the relaxation process via the impurity bound state. 
On the other hand, the peak vanishes at a non-impurity 
site. This is the remarkable difference in T^^ between 
the impurity site and non-impurity site NMR in uncon- 
ventional superconductors. In contrast, the temperature 
dependence of T^^ does not change between the two sites 
for the conventional s-wave. Therefore, observation of 
the peak by the impurity site NMR measurement can 
provide strong evidence of an unconventional supercon- 
ducting state. 
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